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Abstract: One of the most precise measurements of the strong couphng constant Qs(Mz) 
is obtained in the context of global analyses of precision electroweak data. This article 
reviews the sensitivity of different electroweak observables to ag and describes the pertur- 
bative uncertainties related to missing higher orders. The complete renormalisation scale 
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1. Introduction 

High precision measurements of the parameters of the Standard Model (SM) have been 
performed over the last 15 years in particular at LEP, SLC and TEVATRON. Cross sec- 
tions, asymmetries, masses and widths of the electroweak (EW) gauge bosons have been 
determined with a relative accuracy of often better than one per mil. The measurements as 
a whole over-constrain the SM and allow for an internal consistency check. In a global fit 
to these data certain unknown parameters like the mass of Higgs boson can be determined, 
or other quantities not directly measured at LEP like the mass of the top quark can be 
inferred from the LEP data alone [1]. The sensitivity of the electroweak data to these 
parameters arise from higher order corrections. 

The LEP and SLD experiments have carried out global SM fits to combined data 
from various experiments and determine five parameters simultaneously: the masses of the 
Z and Higgs bosons and of the top quark, the hadronic vacuum polarisation Aa^^^ and 
the strong coupling constant as{Mz) constitute one often used set of parameters. The 
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strong coupling constant plays a special role in these fits. It is essentially determined by 
liadronic observables, for which complete next-to-next-to-leading oder (NNLO) calculations 
are available. The measurement of ag benefits from third order 0{a^g) perturbative QCD 
calculation, which is not yet complete for other variables like jet rates, event shapes or 
fragmentation functions. The theoretical systematic uncertainties related to missing higher 
orders are expected to be smaller than for NLO-based determinations, even including all- 
orders resummation used for analyses of event-shape distributions [2]. 

A detailed analysis of the perturbative uncertainty associated to this kind of measure- 
ment of CKs from electroweak observables is the purpose of this article. The renormalisation 
scale dependence of the theoretical predictions is taken as indicative for the systematic the- 
oretical uncertainty originating from missing higher orders, keeping in mind that the real 
uncertainty may in fact be different and will be revealed only once the higher orders are 
actually calculated. Albeit this limitation, the renormalisation scale variation is commonly 
used in other processes and for other variables as well and may therefore at least be used 
for a relative comparison of different measurements. 

Other methods have been proposed in the past to estimate the uncertainty for EW 
observables. Instead of a scale variation certain classes of higher oder terms where cal- 
culated [3] in order to improve the convergence of the perturbation series by minimising 
renormalon effects, and the difference with respect to the standard NNLO result was taken 
as uncertainty. This procedure leads to a very small estimate of ±0.0005 for the uncer- 
tainty of a^Mz, which should be taken with care given that only a subset of the higher 
order terms have been calculated. 

The scale variation method was investigated in [4], where a parameterisation [5] of 
Rz, the ratio of hadronic to leptonic width of the Z boson was used to derive an uncer- 
tainty of +0.0028 —0.0005 for a^. However, the parameterisation embodies only the the 
scale dependence of the massless NNLO part of the calculation but neglects the also scale 
dependent massive quark and mixed EW(8)QCD corrections. 

The objective of this paper is the complete evaluation of the renormalisation scale 
dependence for all relevant observables and an analysis of the contribution from different 
classes of higher order correction to the systematic uncertainty defined by a scale variation. 

The experimental systematic uncertainties for ag are studied in detail by the LEP 
experiments as well as the correlation between and the other EW observables. Taking 
the leptonic cross section at the Z pole as observable from which ag is determined, a 
variation of the Higgs mass in the range from 100 to 1000 GeV entails a change in ag of 
about 0.0022, which has to be compared to an experimental uncertainty of ±0.0030. 

This paper is organised as follows. In Section 2 the EW observables are briefly pre- 
sented, in Section 3 the theoretical predictions for the widths, the effective couplings and 
the final state QCD corrections incorporated in the radiator functions are discussed. In 
Section 4 the theoretical uncertainties of the observables arising from the renormalisation 
scale variation are evaluated. The scale dependence is used in Section 5 to assess the per- 
turbative uncertainty for measurements of Ug extracted from global fits, supplemented by 
an experimental cross-check using test data for the EW observables. The conclusion and 
summary are given in Section 6. 
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2. Electroweaik; observables 



The sensitivity of certain electroweak observables to as arises mainly through pure QCD 
corrections 0{a^) to the decay widths of the Z boson into hadronic final states. In ad- 
dition, mixed QCD^EW corrections 0{aas) to the electroweak couplings give rise to a 
dependence of both hadronic and Icptonic observables on ag. Numerically the former cor- 
rections amount to about three percent, while the mixed corrections, being suppressed by 
a factor of a, are below one per mil. In practice only observables containing the hadronic 
or total widths significantly contribute to the determination of 
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Beyond the width itself the R ratio 
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is of interest as the mixed corrections to the couplings cancel to a large extent in the ratio. 
For measurements of ag the most important observables are the leptonic and hadronic pole 
cross sections 



127r 



r,r 



eJ- h 



M|r| ' 



(2.3) 



In particular exhibits a good sensitivity to ag which allows for a precise measurement 
of Us from this observable alone [1]. As the leptonic pole cross section = a^j Rz is 
not independent of the other variables, it is not included in the global analyses. Realistic 
observables are F^, Rz and (7° in the sense that they are independent observables with a 
substantial sensitivity to a^. 

The calculations presented in this paper are carried out throughout with the elec- 
troweak library ZFITTER version 6.36 [6]. If not stated otherwise, the following numerical 
input values are used: 



Mz = 91.1875 GeV 
mt = 175 GeV 
Mh = 150 GeV 



(2.4) 



l/a(0) = 137.0359895 
a,(Mz) = 0.1185 



3. Pcirtial widths 



The dependence on ag and on the renormalisation scale of the relevant EW observables 
through the widths are given in the following sections. 
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3.1 Dependence of the widths on as 

The partial width of the Z decay to a pair of fermions can be cast into two different 
expressions for leptons and quarks in order to incorporate the different types of radiative 
corrections. The width for lepton pairs / = e, r is given by: 



(3.1) 



while for quark pairs q = u,d, s, c, b another expression is used: 

r, = To Nc \pl\ [bip RliMi) + RliMl)] + A^^^^^^ . (3.2) 
The basic width Fq is the given by 




G 

To = ^ = 82.945(7) MeV , A^c = 3 , (3.3) 
24V27r 

and Nc is a QCD colour factor. In the case of leptons mass effect terms oc are explicitly 
taken into account in Eq. 3.1, for quarks these effects are embodied in the radiator functions 
Ry and i?^, which also account for QCD and QED radiative corrections. The core of the 
sensitivity of the widths and related EW observables to ag stems from the dependence 
of the radiator functions on a^, and their renormalisation scale dependence dominates the 
theoretical uncertainty for a measurement of q.s. The dependence of the widths and derived 
observables is depicted in Fig. 1, where the change of theses quantities normalised to their 
value at as(Mz) = 0.1185 is shown. Since the leading term of the QCD correction is 
1 + as/TT, the dependence of the widths on is basically linear, as higher order terms are 
suppressed by powers of Qs/tt. 

For the leptonic widths a very small change below 0.1 per mil is observed, induced 
by two-loop corrections to the effective EW couplings and to the vacuum polarisation 
contribution to a{s). The hadronic widths exhibit a stronger and opposite dependence 
resulting in a relative change between 4 and 6 per mil over the range of ag between 0.11 and 
0.13, this dependence is induced by the radiator functions. Among the hadronic widths 
there are clear differences between up- and down-type quarks, the d- and the s-quark 
dependences are identical. Furthermore, finite mass quark effects entail small differences 
between the u- and the c-quark, more visible between the d-quark and the b-quark. The 
b-quark behaviour is accidentally very close, but not identical to the sum over all flavours. 

The dependence of the EW observables on shown in Fig. 1 is, as for the widths, 
almost linear and the relative change is between 2 and 8 per mil in the considered range. 
The observable with the strongest dependence is crj^, where the radiator functions enter 
quadratically in the denominator. 

For quarks additional non-factorisable EW^QCD corrections A^^^qj-.^ for the widths 
are not part of the radiator functions. These corrections are numerically very small (less 
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0.11 0.1125 0.115 0.1175 0.12 0.1225 0.125 0.1275 0.13 
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Figure 1: Dependence of the widths (a) and the EW observables (b) on as{Mz)- Shown is 
the normahsed ratio of the observable to its reference value at as{Mz) = 0.1185. The other SM 
parameters are kept fixed to their nominal values. 



than one per mil) and are taken as fixed numbers from [7, 8] 

^Ew/QCD = -0-113 MeV U-, c-quarks (3.4) 
= -0.160 MeV d-, s-quarks 
= -0.040 MeV b-quarks 

The complex- valued variable in Eqs. 3.1, 3.2 measures the overall strength of the neutral 
current interaction in the // channel and the effective coupling g{; can be expressed in terms 
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of the ratio of effective vector and axial-vector couplings 

gl = ^ = l-4\Qf\{4sl + If) , (3.5) 

where defines an effective mixing angle for flavour / with Sw = sin^ 9w given by 

Sw = ^ ~ ^ c^^ = -j^ ■ (3.6) 

The weak isospin /j^^ is ±1/2, the electric charge Qf is +|/ — for up-/down-typc quarks 
and the 0{a^) term ij originating from 77 and Z7 polaristion operators is given by 



r2 „ 2 / ^ 35 



If = a' is) 
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The effective couplings of the Z decay ^-iid p{ incorporate radiative clcctrowcak correc- 
tions up to two loops and their full expressions are given in [6]. The factorisable EW(8)QCD 
corrections 0{aas) shall be studied here, as they induce the ag dependence to the effective 
couplings. 

The running QED coupling denoted by a(s) is given by: 

a(s) = j^. ^ . (3.8) 

The main contribution to the running coupling stems from the hadronic and leptonic 
vacuum polarisation. The leptonic part has been calculated at third order [9] 

Aai^p{M^) = 0.03149767 , (3.9) 

with negligible uncertainties. The hadronic contribution of the five light flavours is related 
via the dispersion relation to Ry (equivalent to Rz in the continuum) from which it can 
be extracted [10] 

A4ad(M|) = 0.02761 ± 0.0036 . (3.10) 

The contribution from the top quark is small but depends on the top quark mass, for 
mt = 175 GeV 

Aa*(M|) = -5.776 10"^ . (3.11) 

An explicit dependence on appears in the 0{aas) correction to a(s) representing gluonic 
insertions in tt loops [11]. Of course the gluon exchange also occurs in light-quark loops 
and is accounted for in the experimental determination of AaJ^^^. The correction for the 
top quark reads 

A«'(MJ) = (Re (i5M) _ 4C(3) + 5) . = ^|±!1 . = 1.2020569 . 

(3.12) 

where the expression for Vi{r) is given in [11]. The numerical value is 

Aa°"-(.s) = -1.02 10"^ . (3.13) 

The dependence of a{s) on ag is very weak, its relative change is about 10~^ for a variation 
of as{Mz) between 0.11 and 0.13. 
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3.2 Effective electroweak couplings 

f f 

The effective electroweak couplings pz and Kz contain various self-energy terms, calculated 

for the EW part at NLO with two-loop corrections and at 0{aas) for the mixed EW^QCD 

corrections, including the terms leading in mt of the ©(aa^) contribution. The expressions 

for EW part are given in [6], here only the terms relevant for the ag dependence are 

f f 

summarised. In a convenient decomposition pz and Kz are split into leading and remainder 
contributions, each being gauge invariant separately. The dominant leading term is re- 
summed to all orders in perturbation theory and the sub-leading remainder is calculated in 
fixed-order theory. The couplings pz and Kz can be expanded in the following combination 
of leading and remainder terms: 
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(3.14) 



(3.15) 



The transformation factor fa accounts for the conversion of couplings a ^ G/j, [12] 
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V2G.Mzs 
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(3.16) 



The three-loop QCD correction to the p parameter arising from top-quark loops is given 
by [13] 



p«- = 3xjc,^ + c,. 



(3.17) 



xt 



Ctl = -2.86, ct2 = -18.18, 



which also includes corrections to the term leading in G^mjaj. derived in [14]. This is the 
dominant QCD correction to pz and Kz and its renormalisation scale dependence determine 
the perturbative uncertainty for pz and k{ derived in Section 4. 

The expansions of p{ and Kz have also the remainders of the renormalisation parameter 
Ar [6] in common, the component containing the QCD corrections is given by 
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ca = 1.409 , ci2 = -12.805 . (3.22) 
The QCD correction for the flavour-dependent remainder of p{ is 

=P^'''' + ^fe^ (3.23) 

I aas{Mz) yj(t) + yj(b) + 2 ^ ^^ ^4) 

yT(g) = 1 - 4Q,s^ (3.25) 
The expression for contains another QCD correction for the remainder 



AC.^^"^ = /^^^^ + TTl^^tW - 3x,c„2 ( ^^"l ' , (3.26) 



Ca2 = 0.644 , (3.27) 

^QCD ^ ^^Mtl ^^^^^ (^^^ ^2) (3_28) 



7r2 2s4 



+ -9 V f^og^VF (3.29) 



The remainder functions drrem, dprem ^'Hd dKrem describing the 0{aas) contribution to 
the bosonic self-energies have been derived analytically in [11]. 

In the case of b-quarks two additional one-loop vertex diagrams, absent for light quarks, 
are generated by the large mass splitting between the t- and the b-quark and contribute to 
the widths [15]. These corrections entail a modification of the Zbb decay amplitude form 
factor [6], which in turn affects the effective couplings. If p' and k' denote the modified 
couplings given in [6] , then the corrected couplings for the b quark are obtained by 

p^z = p' (1 + , (3.30) 

4 = ' (3-31) 



l + Tb ' 



vr , . . rrit 



n = -2xt [1 - -as{mt) + xtT2 j , (3.32) 

where the term in as{mt) was obtained by [16] and the expression for T2 can be found 
in [17]. 

The dependence of the couplings on ag is shown in Fig.2 for the absolute value of pz 
and the squared module of g^, these are the relevant terms for the widths according to 
Eqs. 3.1, 3.2. The dependence of pz on is weak, the relative change in the considered 
range is below 2 • 10^'^. All quark flavours and the leptons exhibit practically the same as 

dependence, except the b-quark for which it is even weaker and opposite. 

The representation of the relative change as ratio is inadequate for since its 

absolute value is rather small, e.g. for leptons jy^P = 0.00555 and jp^l = 1-00517 at 
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Figure 2: Dependence of the \p{,\ (a) and \gl\^ (b) on as(Mz). For \pl\ the normalised ratio and 
for l^zp the shift of the coupHng to its reference value at as(Mz) = 0.1185 is shown. 



as{Mz) = 0.1185. Instead the shift defined by \gi\^{as) - \gi\'^{as = 0.1185) is a better 
indicator for the dependence on ag. The absolute change of l^^p is between 1.4 • 10~^ for 
up-typc quarks and 0.4-10"^ for Icptons. This translates into a relative change much larger 
for leptons than for quarks, given the small size of Igz]"^ for the leptons. 

The observed as dependence of the leptonic widths shown in Fig. 1 is dominated by 

f 

dependence of /O^ on a^. In the case of quarks, however, the couplings contribution to the 
as dependence of the widths is sub-leading, the main properties are determined by the 
QCD final state corrections in the radiator functions. 
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3.3 Radiator functions 



Final state QCD and QED vector and axial vector corrections to the quaxkonic widths 
Eq. 3.2 are embodied in the radiator functions 

^ - 1 ^ 3 2a(s) , «.(g) _ 1^2"(-'^) (^sjs) 



vr TT 



+ 



Co2 + Co ( — o 



m; 



TT 



TT 



, rnl{s)+ml{s) ( a^s) 

H <-y23 



TT 



s 

-.4/ 



-•21" 



TT 



TT 



23 



TT 



C42 - In 



m2(s) 



2 1 ^41 

TT 



TT / 

cr2+cr2^in^ 



C42 - In 



m2(s) 



TT 



TT 



+12- 



TT 



16 



155 + 6 In 



m2(s) 



TT 



(3.33) 



4^7r TT 4^7r vr 



+ 
+ 



Cos 

2/' 



(2/f>) 



T 



(3) 



m: 



t / J 



TT 



i^f" 



23 



TT 



+ 



20 T" "-^21 Z r '-^22 



TT 



TT 



mo(s) 



8 1 , mf 

1 In — 

81 54 s 



IT 



C42 — In 



TT 



C42 — In 



TT 



+ 



S ^40 "T '^41 " 



+ 



TT 



C4l + C//In^ 



vr 



TT 



(3.34) 



Finite mass corrections are retained only for the b- and c-quark, i.e. rUg = for q=u,d,s, 
and the terms mq{s) represent the running quark masses in the MS scheme. The term 
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denotes the other quark mass in doublet, it is nic for q=b and mh for q=c. The different 
terms of Eq. 3.33 and Eq. 3.34 and their coefficients can be organised in the following 
classes of corrections. 
Massless non-singlet corrections [18-21]: 
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with the number of active flavours n/. 
Quadratic massive corrections [22]: 
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Quartic massive corrections: 
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Power suppressed t-quark mass correction: 
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Singlet axial corrections: 
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Here, the Riemann Zeta function C, is defined by 
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with particular values 

C(2) = 1.6449341 , C(3) = 1.2020569 , C(5) = 1.0369278 
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The evolution of the radiator functions with Ug is shown in Fig. 3. The vector radiator 
functions for the different flavours are very similar and can barely be distinguished. Their 
dependence on is almost linear and they are increasing by 0.6% for a change of q,, from 
0.11 to 0.13. The axial-vector radiator function exhibits a prominent flavour dependence. 
For up-type quarks R\ increases by 0.8%, for down-type quarks the change of the axial- 



vector radiator function is only of 0.4%. The third component of the weak isospin /, 
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Figure 3: Dependence of the radiator functions (a) and (b) on as{Mz). The ratio of the 
radiator function at a given value of as{Mz) and the reference value at as{Mz) = 0.1185 is shown. 



present in the singlet axial component of R\, generates the up/down-quark difference. On 
top of this difference the b-quark mass corrections are important and entail a change of R\ 
of 0.6%. 

3.3.1 Running quark masses 

The running quark masses mq{s) in the MS scheme are related to the fixed- valued pole 
masses Mq. For the c-quark at s = M| and = 4 it follows 
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The running mass is evolved according to the renormahsation group equation from the scale 
of the pole mass to the scale of the process s in a two-step evolution M| — > — >■ s: 
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For the running b-quark mass the same procedure is applied with a single evolution from 
M? to s. The coefficients in Eq. 3.62 are given by: 
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The coefficients of the Beta and Gamma functions are: 
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(3.68) 
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The renormalisation scale dependence of the couphng constant can be parameterised at 
3-loop level as function of 
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Technically, for a given input value of as{Mz), Eq. 3.71 is solved numerically for 
n/ = 5 in order to obtain A^. The scale parameters for n/ = 4 and n/ = 3, A^ and 

Alfi., required for the evolution of ag to the scales of the quark pole masses Mi, and Mc, 
are derived using the matching condition 
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3.4 Quantifying higher order contributions 



(3.72) 



(3.73) 



The theoretical prediction for the widths and derived EW observables consists of three basic 
ingredients, each incorporating different classes of higher order contributions: the NNLO 
massless terms, the quark mass corrections and the mixed EW(8)QCD corrections. For the 

evaluation of associated perturbative uncertainties it is essential to study the impact and 
size of the these contributions to the full theory. The sensitivity of widths and realistic 
observables to the three classes of corrections is illustrated in Table 1, where the relative 
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Theory 



u 

To 
Rz 



NNLO 



300.057 
382.901 
299.994 
375.807 
1741.66 
83.796 
2495.08 
41.4798 
2.0002 
20.737 



MeV 
MeV 
MeV 
MeV 
MeV 
MeV 
MeV 
nb 
nb 



nominal value 



NLO without quark without mixed 
mass corrections corrections 



-0.13 
1.28 

-0.18 
1.23 
0.77 
0.00 
0.54 

-0.31 

-1.08 
0.78 



0.00 
0.00 
0.21 
4.13 
0.93 
0.00 
0.65 
-0.37 
-1.29 
0.93 



relative change [% 



'00 



1.66 
1.45 
1.67 
1.47 
1.53 
1.19 
1.37 
-0.03 
-0.37 
0.34 



Table 1: Nominal vahics of widths and EW obscrvablcs and their relative change in per mil 
observed when certain classes of corrections are omitted. The calculations are carried out for 
Q!^(Mz) = 0.1185. 



change of the full prediction with respect to downgraded calculations neglecting certain 
terms is given. The change from NNLO to NLO for the quarkonic widths is between one 
per mil for down-type quarks and 0.1 per mil for up-type quarks, this large difference is 
generated by flavour-dependent contributions to R\ and (/|. Neither NNLO nor quark 
mass corrections have any sizeable impact on the leptonic width, but the mixed EW^QCD 
corrections entail changes of 1.5 per mil to all widths, larger than the final state NNLO 
corrections. Not surprisingly, the quark mass corrections are as large as 4 per mil for 
b-quarks but drop to 0.2 per mil for c-quarks. The total width absorbs an average effect 
induced by the widths and the other realistic observables are in general less sensitive to 
these corrections cancelling in the ratio of widths. Among the realistic observables cj^ is 
least and most sensitive to the higher order corrections, their size ranges between 0.3 
and 1.3%o. 



4. Theoretical uncertainties for EW observables 

The sensitivity of a given electroweak observable to originate on one side from the 
QCD corrections incorporated in the radiator functions and on the other side from the 
mixed EW(X'QCD corrections to the effective couplings. A measurement of Ug using EW 
observables is subject to a systematic uncertainty stemming from missing higher orders 
in the perturbation series. The yet uncalculated higher orders are inherently difficult to 
access. A conventional method of estimating the perturbative uncertainty consists of a 
variation of the renormalisation scale /i. The natural scale of the process is usTially taken 
to be y/s in c"'"c~ annihilation, and subsequently in the case of Z peak obscrvablcs ji is set 
to Mz- Neither this particular choice for the nominal scale nor the range of variation for 
H are unambiguous [2]. Following the convention applied in analyses of e+e~ event-shape 
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variables, the perturbative uncertainty is estimated by changing = lil\fs in the range 
1/2 < < 2. 

A variation of the renormahsation scale induces a change of the value of ag as given 
in Eq. 3.71. At NLO and beyond this change is compensated by a modification of the 
(N)NLO terms, resulting in a residual dependence at (N)NNLO. The details of the scale 
dependence of the radiator functions is discussed below. 

The mixed 0{a.as) corrections are complete only at leading order in a^, but the dom- 
inant three-loop correction in Eq. 3.17 leading in is included. For these O(aa^) terms 
the explicit scale dependence is taken into account. 

Given the overall small size of the EW(8)QCD corrections, they do not contribute 
significantly to the scale dependence of the realistic observables. 

4.1 Renormalisation scale dependence 

Dimensional regularisation introduces a renormalisation scale ^ at which the coupling con- 
stant is defined. Thereby, the coefficients in the expansion of -R^y^ acquire an explicit 
dependence on this scale, which is only at all orders completely compensated by the scale 
dependence of as{\i)- For a NNLO calculation, the residual scale dependence is N^LO. The 
nominal value of the /x scale is set to the scale of the process = s. 

The expression for Ry^j^ in Eqs. 3.33, 3.34 are valid only for /i^ = s, for different 
renormalisation scales terms proportional to powers of In/Lt^/s appear. For a generic power 
series of the type 

n 

H=i:^(^)", (4.1) 

1=0 

the NLO coefficient C2 becomes a function of fj, 

C2 ^ C2(At) = C2 + /3oCi In ^— . (4.2) 

s 

It is important to note that two quantities depend on the renormalisation scale in Eqs. 3.33 
and 3.34: the coupling constant as{^) (Eq. 3.71) and the running masses mq{p,) (Eq. 3.62). 

In order to to simplify the formulae for the scale dependence of the radiator functions, 
it is convenient to re-order Eqs. 3.33 and 3.34 in terms of powers of the running masses: 

f=<l,g',c,b 1=0 j=Q 

where a. = For each quark not only the mass of the actual quark, but also the 

mass the b- and c-quark masses intervene in the radiator functions, each with different 
coefficients. In the MS scheme the expansion of the scale evolution of powers of the running 
masses reads as 



ml{s) = m2(/i) 1 + 27oLa, + {joPqL'^ + 2jiL + 2^1l^) 



(4.4) 



a. 



ml{s) = mq{ii) (1 + 67oLas) , 



(4.5) 
(4.6) 



with L = Inx^. The new coefficients df'^ are related to those of Eqs. 3.33 and 3.34 by the 
following formulae. 

Massless terms: 
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(4.7) 
(4.8) 

(4.9) 

4.10) 
4.11) 
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4.13) 
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,d 



16,? 



(4.25) 
(4.26) 



Finally, to get the renormalisation scale dependence of the radiator functions, the terms of 
the type mq{s)'^-dijal{s)/TT in Eq. 4.3, dropping for clarity the axial- vector/vector and 
flavour indices, have to be replaced by the following expressions: 



j:d 

j=0 



oj- 



TT 



doo + dio^s + (c?20 + dioPoL) ag 
+ [dso + {dioPi + 2d2o/?o) + dioPW) o^s^, 



m, 



\{s)Y,d 

3=0 



ai{s) 



TT 



do2 + {di2 + 2do2loL) 



(4.27) 



(4.28) 



+ {di2f3oL + d22 + 2di2joL + do2lof3oL + 2do2liL + 2do2loL 
+ {^i20lL^ + di2/3iL + d32 + 2d22poL + ^do2loPiL^ + do27oPiL'^ 

+2do2lifioL^ + 2do2l2L + 2do2llPoL^ + ^do2loliL^ + ^do2llL^ 



+3di2loPoL^ + 2d22loL + 2di2liL + 2di2llL^ « 



,j4- 



c4(s) 



3=0 



TT 



doi + {dii + 4(io47o^) "s 



(4.29) 



{di4^f3oL + ^24 + 2doAloPoL^ + 4cio47i-^ + 8do47o-^^ + 4di47oL) 



m. 



':{s)Y.d,,^^mii„) 

j=0 



TT 



do6 + {die + QdoejoL) 



(4.30) 



The dependence of the radiator functions on the logarithm of the renormalisation scale 
Inx^ for a fixed input value of as(Mz) is shown for each quark flavour in Fig. 4. 

The shape of the scale dependence of Ry is almost identical for all flavours, except 
a small quark mass modification for the b-quark. A maximum in Ry appears around 

lnXfj_ = 0.5 and a minimum at —1.75, spanning a difference of one per mil. The overall 
scale dependence of the axial- vector component i?^ is twice as large as the one of Ry. As 
already observed for the dependence on Ug in Fig. 3, the shape of R^ is clearly different for 
up- and down-type quarks. For up-type quarks i?^ becomes maximal at Inx^ = —0.3, the 
maximum is close to the nominal value at Inx^^ = 0. The shape of the scale dependence 
is opposite for down- type quarks: a minimum appears at Inx^^ = —1.1. Considering the 
range of variation for from 1/2 to 2, corresponding to a range for Inx^ from —0.7 to 
0.7, it appears that largest deviation from the nominal point at Inx^ = is not always 
obtained at the endpoints, but sometimes at smaller variations. Therefore, when assessing 
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Figure 4: Dependence of the radiator functions on the renormalisation scale. The vector correction 
Rv (a) and the axial vector correction R'i (b) are shown for the quarks q = u,d, s, c, b as function 
of InXfj,, normalised to their value at a;^ = 1. 



the uncertainty for the observables studied in the following, the endpoints of Inx^ (lna;~ 
and lnx+) have been chosen to correspond to the largest change in the observables, within 
the pre-defined range | Inx^l < 0.7. 

The dependence of the effective couplings on Inx^ is shown in Fig. 5. The quantities 
\pz\ and l^zp, as they appear in the expressions for the widths, are shown as ratio for 
but as shift for g^, i.e. |5zP(ln.T^) — \gz\'^{lnXf^ = 0). The shape of the couplings scale 
dependence is rather different from the radiator functions dependence and less structured, 
given by the interference of complete two-loop and incomplete leading three-loop corrections 
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Figure 5: Dependence of the effective couplings \pz\ (a) and \gz\'^ (b) on the renormaHsation scale. 



f 

to the mixed eorrections. For all light quarks and leptons the relative change of |/3^] is only 
about 10^^, for the b-quark mueh larger and amounts to 0.8 per mil. The absolute change 
of |(7^P is about 2 • 10~^ for light quarks and 0.6 • 10~^ for leptons and slightly less for 
b-quarks, corresponding to a relative change of less than two per mil for all quarks but 
almost one percent for the leptons. 

Turning to the widths, for hadronic final states both the effective couplings and the 
radiator functions discussed above contribute to scale dependence. For the Icptonic widths 
only the effective couplings depend on the renormalisation scale through the mixed cor- 
rections. The relative magnitude of contributions is given by the formulae for the widths, 
Eq. 3.1 for the leptons and Eq. 3.2 for the quarks. The scale dependence of the widths and 
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Figure 6: Dependence of the widths (a) and selected realistic observables (b) on the renormalisation 
scale for a^Mz) = 0.1185. 

of the realistic observables are shown in Fig. 6. 

The shape of the scale dependence for the widths depends on the fermion type. For 
the u- and c-quark a maximum is found close to Inx^ = 0, and variations of Inx^ in any 
direction entail a decrease of the width. The partial widths into d- and s-quarks increase 
monotonically from Inx/j, = —0.7 to Inx^ = 0.7, similarly for the width of the b-quark, 
albeit with a flatter shape. The total hadronic widths of the Z boson emerges as sum of 
the quarkonic contributions, leading to an average shape of its scale dependence with a 
minimum at —1.6, a maximum at 0.7 with a difference of one per mil between them. 

The scale dependence of the leptonic widths is clearly much weaker than that of the 
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quarkonic counterparts. In the central range of In the leptonic width changes by about 
a tenth of the hadronic width's change. The sensitivity of to the renormahsation scale 
arise through the EWxQCD corrections, essentially those incorporated in p^. 

4.2 Perturbative uncertainties for observables 

The perturbative uncertainty for widths and realistic observables are defined as the dif- 
ference between the observables value at Xfi = 1 and at , for a given input value of 
as(Mz). The range of variation for z^j, defined by and x~, is chosen to generate the max- 
imum and minimum values of the observables within the pre-defined range 1/2 < x^ < 2. 
The absolute size of the perturbative uncertainty depends on the input value of a^, for an 
observable calculated at NNLO its uncertainty is scaling with a^. The relative systematic 
uncertainty for the widths and observables is shown in Fig. 7 as function of as{Mz). Here 
and in Table 2 the relative uncertainties for a generic observable O are defined by 



The positive and negative systematic uncertainties are generally asymmetric, for the 
widths of the u- and c-quark only the negative uncertainty contributes, given the particular 
shape of their scale dependence. The uncertainty for the width of down-type quarks is 
typically ±0.05%, for up-type quarks it is —0.03% and for the leptonic width a factor of 10 
smaller ±0.005%. Among the realistic observables cr° has the largest uncertainty of about 
+0.05%, the other variables uncertainty is about half that size. 

For selected values of as{Mz) the uncertainties for the observables are given in Table 2, 
where the scale variation endpoints listed in the last column. 

4.3 Contributions to the scale dependence 

The three classes of higher order corrections analysed in Section 3.4 contribute to the 
perturbative uncertainty of the observables via their different evolution under the renor- 
mahsation scale variation. The size of the NNLO, quark mass and mixed EW(8)QCD 
corrections, given in Table 1, are typically at the level of one per mil, while the perturba- 
tive uncertainties are about 0.1 per mil (see Table 2). The absolute size of the corrections 
alone is not a reliable indicative of their contribution to the uncertainty, which must be 
evaluated from their scale dependence. 

The scale uncertainty is therfore re-calculated without the NNLO contribution (i.e. 
dropping terms in 0{ag) in Eqs. 3.33, 3.34), in which case also the evolution of as(/x) in 
Eq. 3.71 has to be performed at NLO. 

The contribution from the quark mass corrections is estimated by switching off the 
explicit scale dependence of the running quark masses (i.e. setting mq{s) = mq{fi) in 
Eqs. 4.5 - 4.6). 

The impact of scale dependence of the mixed EW(X'QCD corrections is tested by elim- 
inating the terms in the 0{aa'l) corrections depending explicitly on X/^. 



AO = 




(4.31) 



0{x, 



1) 
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The scale dependence obtained under these conditions is compared to the nominal 
scale dependence of the complete NNLO prediction in Fig. 8, considering as example Tz 
and cr^. 

The change from NNLO to NLO entails a strongly increased scale dependence for the 
quarkonic widths, but with opposite effect on up- and down- type quarks, resulting in a 
weaker enhancement for the realistic observables. The uncertainties for the observables 
obtained as before by a scale variation between and x~ are summarised in Table 3 for 
the theories with modified scale dependence of higher order corrections. For the realistic 
observables the downgrade from NNLO to NLO results in an increase by a factor of two of 
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Table 2: Systematic perturbative uncertainties in per mil of EW observables for different values 
of as(Mz). The last colum indicates the values of the renormalisation scale a;^ corresponding to 
the maximum variation of the observables with respect to their nominal values at a;^ = 1, within 
the pre-defined variation range 1/2 < a;^ < 2 

the perturbative uncertainty and the asymmetry between positive and negative uncertainty 
is also enlarged. For realistic observables the shape of the scale dependence without the 

running masses is similar but steeper than in the NNLO case. It has a large impact only 
for the b-quark, where the uncertainty is increased by a factor of three, which propagates 
to an increase of about 50% for the realistic observables. 

The effect of dropping the scale dependence of the mixed corrections has in contrast 
to the other variations a net impact on the leptonic width, since this is the only source 
of sensitivity to QCD effects. The perturbative uncertainty for Vi is increased by a factor 
of three, the quarkonic width uncertainty by about 30-40 %. The realistic observables of 
ratios of widths are less sensitive to the mixed corrections and a moderate enhancement 
of the uncertainty of a few percent is observed, while for Vz the effect from quarkonic and 
leptonic widths adds up to an almost doubled uncertainty. 

In conclusion the stability of the predictions under scale variations depends crucially 
on the NNLO corrections, but depending on the observable also quark mass and to lesser 
extent mixed corrections contribute significantly to the accuracy of the calculations. 
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Figure 8: Renormalisation scale dependence of Vz (a) and a\ (b) using the full NNLO prediction 
compared to reduced theories at NLO only, neglecting the scale dependence of the running quark 
masses and of the mixed EW(8)QCD corrections, for as(Mz) = 0.1185. 



5. Pert ur bat ive uncertainties for 



In the context of global analyses of world electroweak data [1] q;s(Mz) is fit together with 
four other free parameters of the standard model: Mh, Mz, mt and Aq:|^^^^. The correlation 
between OLg and the other parameters is small. The obscrvablcs included in the fit with a 
sizeable sensitivity to ct^ arc i?^, V z and a^. The information from the leptonic pole cross 
section is included in the other observables. This particular selection of observables is 
the result of an optimisation for the best accuracy for the five free SM parameters, which 
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Table 3: Perturbative uncertainties (in per mil) for selected observables obtained from renormali- 

sation scale variation at as(Mz) = 0.1185. The nominal NNLO uncertainty is compared to reduced 
theories where the scale dependence for certain classes of corrections is switched off. 

may not necessarily be the optimal for ag alone when the other parameters are fixed to 
their SM values. For example has actually the best sensitivity through the inverse 
squared radiator functions and may be used alone to determine from a single parameter 
fit, avoiding thereby the otherwise required correlations between the observables, and the 
dependence of ag on the Higgs or top mass may be investigated. 

Having in mind this scenario, the procedure is to determine in a first step the pertur- 
bative uncertainty for measurements of ct^ using single selected observables and to estimate 
in a second step the uncertainty for a global fit including several variables. 

The basic principle for the uncertainty estimation was developed in [2]: the systematic 
uncertainty for a given observable and fixed value of as{Mz) (i.e. as obtained from a 
fit to this observable) is evaluated in the theory by variation of the renormalisation scale 
x~ < X/j, < x^. The change of the observables under the scale variation can also be 
generated by a variation of the input value of at fixed a;^ = 1. This procedure leads 
in general to two alternative values of ag corresponding to the changes of the observable 
for xTT and xt- The difference between the nominal value of a, and these two alternatives 
finally determine the perturbative uncertainty of a^. 
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Figure 9: Positive and negative contributions to the perturbative systematic uncertainty oias{Mz) 
determined from various partial widths of the Z boson. The uncertainty using Ti obtained at NLO 
is divided for representation by a factor of two. 

The uncertainty itself depends on the type of observable and on the input value of 
ag. At NNLO the size of the perturbative uncertainty scales as . It is instructive to 
consider the uncertainty for for the pseudo-observable widths in a first step, in order 
to understand their contribution to realistic observables. The systematic uncertainties are 
shown in Fig. 9 as function of the input value of in a relevant range from 0.11 to 0.13. 

As expected from the scale uncertainty of the width itself, there are large differences 
between the uncertainties of Og determined using the widths of up-, down- type quarks and 
leptons. The size of the uncertainty for ag is between one and two percent for quarks and 
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about 4 percent for leptons, QCD corrections for the leptonic widths being calculated only 
at NLO. Given the shape of the scale uncertainty for the width of the u- and c-quark, 
the resulting uncertainty of is essentially one-sided. Also for the other widths a certain 
asymmetry in the uncertainty is observed, the positive (upward) uncertainty is generally 
larger than the negative (downward). This asymmetry may well be a technical artefact of 
the scale variation prescription, and conservatively the maximum of the positive and nega- 
tive uncertainty is assigned as a symmetric uncertainty. For selected input values of as{Mz) 
the symmetrised uncertainties are given in Table 4. The uncertainty of as determined from 



as{Mz) 


0.110 


0.115 


0.120 


0.125 


0.130 


Tu 


0.00049 


0.00058 


0.00068 


0.00079 


0.00092 




0.00154 


0.00189 


0.00229 


0.00277 


0.00333 


Tc 


0.00047 


0.00056 


0.00066 


0.00077 


0.00089 




0.00066 


0.00085 


0.00109 


0.00136 


0.00170 




0.00083 


0.00102 


0.00123 


0.00147 


0.00176 


ri 


0.00323 


0.00363 


0.00407 


0.00453 


0.00502 



Table 4: Systematic perturbative uncertainties for measurements of from EW observables. The 
symmetric uncertainty is given by the maximum of upward and downward uncertainties obtained 
by a renormalisation scale variation. 



realistic electroweak observables is shown in Fig. 10. The perturbative uncertainty ranges 
between 1.3% for Tz and 1.0% for a^. The largest uncertainty is observed in the case of 
Tz, which is directly proportional to the product of effective couplings and radiator func- 
tions. The uncertainties from the other variables are very similar, between 1.0% and 1.1%. 
In Rz the widths appear linearly in nominator and denominator, in af and quadratic 
combinations of leptonic, hadronic and total widths of the Z boson interplay. As a conse- 
quence the scale dependence of the effective couplings and/or the radiator functions cancel 
to some extent in the ratio. 

Several observables are included in the global EW fits [1], but only Rz, Tz and a^^ 
have a sizeable sensitivity to a^. Effectively, these three observables determine ag and the 
perturbative uncertainty of ag is bound to be an average of their individual contributions. 
The exact weights of each variable to the determination of ag and subsequently to its 
perturbative uncertainties can not be determined in the framework of the present work. In 
order to derive nonetheless an estimate of the perturbative uncertainty from a global fit, a 
new variable 

A = u;iTz + uJ2(Th + ^3Rz , (5.1) 

is introduced in order to approximate the sensitivity of the three real observables to ag in 
a global fit by a single observable. The perturbative uncertainty for ag obtained from the 
A variable is an indicative for the true uncertainty from a combined fit. The weights uii 
are determined from the sensitivity 
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Figure 10: Systematic perturbative uncertainties for measurements of using the EW observables 
Tz, Rz, cr° and as function of the input value of as{Mz). The quantity A denotes a weighted 
average of the uncertainty for Tz, Rz and aJJ, relevant for global analyses of electroweak data. 



where Oj are the three observables, cr^as) is the uncertainty on as from the fit, taken to 
be 0.003 for all observables, and a{0^^^^) is the measured experimental uncertainty for the 
observables themselves. This procedure yields uji = 0.36, L02 = 0.23 and L02 = 0.41. The 
resulting uncertainty from A is also shown in Fig. 10. The upper and lower bounds of the 
combined uncertainty can be estimated assuming that the true weights of the observables 
in a global fit are all positive. In this case the combined uncertainty can not be larger or 
smaller than any of the individual contributions from the single observables. 

In Table 5 the symmetric perturbative uncertainties are summarised for the realistic 
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observables and for the A variable in a narrow range around 0.119. The possible variations 
for the combined uncertainty are also given as A"*" for the upper and A~ for the lower 
bound. As a stability test, the weights in the calculation of A have been set constant for 
all observables. The resulting combined uncertainty does not change by more than 0.00001. 





0.118 


0.119 


0.120 
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0.122 
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Rz 


0.00119 
0.00095 
0.00098 

0.00099 


0.00124 
0.00098 
0.00101 
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0.00102 
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0.00110 
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A 

A+ 

A- 


0.00105 
0.00119 
0.00095 


0.00109 
0.00124 
0.00098 


0.00114 
0.00128 
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0.00118 
0.00133 
0.00106 


0.00122 
0.00138 
0.00115 


-1.50 • 10"^ 


6.22 



Table 5: Systematic perturbative uncertainties for measurements of from EW observables. The 
symmetric uncertainty is given by the maximum of upward and downward uncertainties obtained 
by a renormalisation scale variation. The variable A represents the combined uncertainty for a 
global fit including Vz, o\ and o-°, with upper bound A+ and lower bound A~. 

For a value of agiM'i) = 0.12 the perturbative uncertainty ranges from ±0.00102 for 
g\ to ±0.00128 for Tz, with a weighted average of ±0.00115 (unweighted ±0.00114). 

The dependence of the symmetric uncertainty on the input value of a^iMj) can 
smoothly be parameteriscd to form a + 6 • as(-^z)^- The parameterisation, valid for 
0.11 < as(Mz) < 0.13, allows for an interpolation of the uncertainty between calculated 
points. The parameters a and h are given for each observable in Table 5. The result of 
the parameterisation is compared in Fig. 11 to the exact calculation, which is reproduced 
to good accuracy. 

5.1 Experimental tests 

In order to substantiate the uncertainty estimate, a fit is performed to a set of test data 
taken from the preliminary LEP combination [1]. In the fit to the test data, the value of 
Us is determined for x^^ = 1 (nominal case) and x'l^,x~ for the uncertainty. Only realistic 
observables are considered and the other SM parameters are fixed to the values given 
in Section 2. The central values, experimental systematic uncertainties and correlations 
between the observables used are given in Table 6 (taken from [1]). A fit to the test 



observable 


central experimental 
value uncertainty 


correlations 
rz ctI Rz 


rz[GeV] 

Rz 
af [nh] 


2.4952 0.0043 
41.540 0.037 
20.767 0.024 
2.0003 0.0027 


1.0 -0.297 0.004 
1.0 0.183 
1.0 



Table 6: Test data for the EW observables. The numerical values are preliminary results taken 
from [1], since is not included in the global fit its correlation data is not available. 
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Figure 11: Symmetrised systematic perturbative uncertainties for measurements of Ug using the 
EW observables. The quantity A denotes an weighted average of the uncertainty for Tz, Rz and 
cr^, relevant for global analyses of electroweak data. The symbols show the exact calculation and 
the lines represent the result of the parameterisation. 



data yields results for ag summarised in Table 7. The central values for fits using single 
observables are close to the world average [4] apart from which yields a very low value 
of as with a large experimental uncertainty of 6%. As expected, the smallest error of 3% 
is obtained for af having the best sensitivity. The perturbative uncertainties obtained in 
this test are in good agreement with the data-independent method summarised in Table 5. 
The case of cj^ yielding a small theoretical uncertainty originating from its small central 
result illustrates the bias appearing when the uncertainty is determined from an actual 
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measurement. This bias can be reduced by taking a combined world average for in the 
theoretical estimation of the perturbative uncertainty. 

The experimental precision of fits to several observables is increased to 3% for the 
canonical ensemble of Tz, cr|J and Rz which is also used by the global EW fits [1]. The 
perturbative uncertainty is again consistent with the result obtained previously for the A 
variable. The observables crJJ with less sensitivity may be exchanged by af in the data 
sample to be fit. The correlation of experimental uncertainties between af and the other 
variables is unknown, as a crude approximation the same correlation as for af^ is assumed. 
Under this reserve the experimental uncertainty is further reduced to 2% with essentially 
unchanged perturbative uncertainties. As the experimental uncertainty dominates the 
precision for ag, the replacement of (T° by af might be considered for the global EW fits. 



included 


central 


experimental 


•^n 




perturbative 


observable 


value 


uncertainty 






uncertainty 


Tz 


0.1174 


0.0041 


0.1170 


0.1186 


-0.00043 + 0.00122 


-1 


0.1076 


0.0065 


0.1074 


0.1082 


-0.00015 + 0.00063 


Rz 


0.1231 


0.0037 


0.1228 


0.1244 


-0.00033 + 0.00125 




0.1187 


0.0030 


0.1184 


0.1197 


-0.00027 + 0.00104 


Tz, ctO, Rz 


0.1191 


0.0027 


0.1188 


0.1203 


-0.00043 + 0.00113 


Tz, af, Rz y 


0.1202 


0.0021 


0.1199 


0.1213 


-0.00031 +0.00114 



Table 7: Fit results for as{Mz) from fits to various test data sets. The central result is obtained 
for Xfj, = 1, the perturbative uncertainty from a variation for the renormalisation scale from a;+ to 
x~. )* In the combined fit to Tz, erf and Rz the correlation for af is assumed to be the same as 
for al 

The LEP electroweak working group has adopted another strategy to incorporate the 
QCD uncertainties in the global fits. Prom the perturbative uncertainties for the observ- 
ables Tz, af^ and Rz themselves, given in Table 2, a covariancc matrix was constructed and 
added to the other covariance matrices related to statistical and further systematic uncer- 
tainties. The total covariancc matrix is then included in the global fit. While the central 
value of Os does not change significantly when the QCD covariance matrix is added, its 
total uncertainty does increase and by quadratic subtraction of the non-QCD uncertainties 
a perturbative uncertainty of ±0.0010 is obtained [23]. This result confirms the estimate 
of perturbative uncertainty based on theoretical considerations for the A variable given in 
Table 5. 

6. Conclusions 

A new method has been presented for the perturbative uncertainties at NNLO of measure- 
ments of as obtained from global analyses of precision electroweak data. The systematic 
uncertainties are obtained by a variation of the renormalisation scale in the calculations of 
final state QCD and mixed QCD(giEW corrections for the electroweak observables included 
in the global analyses used to determine ag and other Standard Model parameters. The 
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NNLO massless corrections, quark mass corrections and mixed QCDi^EW corrections are 
observed to contribute to the theoretical predictions for the observables at the level of one 
per mil. Individual contributions to the renormalisation scale dependence have been stud- 
ied in detail and the resulting uncertainty has been calculated for the widths of the Z boson 
into different quark flavours and leptons. For electroweak observables used to constrain 
the Standard Model and determine as(Mz) the corresponding perturbative uncertainty for 
a value of as(Mz) = 0.12 is estimated to be between ±0.0010 and ±0.0013, in average 
±0.0011. The size of the perturbative uncertainty has been cross-checked with a toy fit to 
experimental test data. The determination of at NNLO from electroweak data is one of 
most precise measurements of as{Mz). Its precision is dominated by experimental effects 
yielding a relative uncertainty of about 3%, while perturbative uncertainties contribute 
only 1%. 
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